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INTRODUCTION 

One area of research which offers a great deal of hope that an inexpensive, systematic, 
control system synthesis procedure applicable to a large class of problems can be 
developed is the application of optimal control theory to time-invariant linear systems 
with quadratic performance criteria. For this class of problems, a library of subrou- 
tines has been developed at Ames Research Center^ 1 ) which can quickly perform all of 
the required computations. 

The synthesis procedure presented in this paper is based on the solution of the output 
regulator problem of linear optimal control theory for time-invariant systems. ( ) By 
this technique, solution of the matrix Riccati equation leads to a constant linear 
feedback control law for an output regulator which will maintain a plant in a particu- 
lar equilibrium condition in the presence of impulse disturbances. 

The use of this technique as a basis for flight-control system design has many very 
appealing advantages over root locus and frequency domain techniques. Some of these 
advantages are that this technique is completely automated after selection of a few 
parameters in a performance criterion; can easily handle multiple input-output systems; 
and is generally applicable to any controllable and observable system whose equations 
can be nut into linear, constant coefficient state space form. Further advantages are 
that the resulting feedback control law guarantees system stability ; \ 3 ' requires no 
dynamic compensation; and is optimal in the sense that it minimizes a quadratic cost 
function. 

On the other hand, some fundamental difficulties must be overcome. One of these is 
that the regulator tends to minimize the effects of any disturbance, including desired 
state changes attempted by the pilot. To avoid the necessity of disengaging the regu- 
lator while performing a maneuver, thereby losing the advantage of the augmented sta- 
bility, a feedforward control law must be computed which can be superimposed upon the 
regulator and which will result in a desired maneuver on command. A second difficulty 
is that the control law resulting from the solution of the matrix Riccati equation 
requires measurement of the entire state vector. If this is not possible or desirable, 
then either the approach must be abandoned or an adequate estimate of the unmeasurable 
states must be constructed. This can be accomplished satisfactorily, and with a mini- 
mum order dynamic compensator, with a Luenberger observer. ( ** » However, construction 
of an observer for high order systems can be a very tedious exercise that has, in the 
past, been carried out largely by hand. Before observer theory can be included in a 
general systematic synthesis procedure, an algorithm must be found to construct 
observers that can be carried out in a straightforward manner by a digital computer. 

This paper presents two simple algorithms that can be used in an automatic synthesis 
procedure for the design of maneuverable output regulators requiring only selected 
state variables for feedback. 

The first algorithm is for the construction of optimal feedforward control laws that 
can be superimposed upon a Kalman output regulator and that will drive the output of a 
plant to a desired constant value on command. 

The second algorithm is for the construction of optimal Luenberger observers that can 
be used to obtain feedback control laws for the output regulator requiring measurement 
of only part of the state vector. This algorithm constructs observers which have mini- 
mum response time under the constraint that the magnitude of the gains in the observer 
filter be less than some arbitrary limit. 
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The algorithms can he useful for the design of flight controllers that are required 
to track a step command input using only selected state variables for feedback, but 
where zero steady-state error (i.e. , integral feedback) is not a requirement. In many 
practical cases , the pilot can be expected to perform this function. 

Illustrative examples of the use of both algorithms for the design of flight controllers 
for the Augmentor Wing Jet STOL Research Aircraft are presented. 

1. SOLUTION OF THE CONSTANT OUTPUT TRACKING PROBLEM 


Let a controllable and observable plant be represented by 

x(t) = Fx(t ) + Gu(t) ; x(t Q ) = x Q (l.l) 

y(t) = Hx(t) (l.2) 

where x is an n-dimensional column vector defined as the state; y, a p-dimensional 
column vector defined as the output; and u, an m-dimensional column vector defined 
as the control or input to the plant. The matrices .F, G, and H are constant, rec- 
tangular, and have appropriate dimension. 

We would like to derive a control law with the properties of an optimal output regula- 
tor and, in addition, with the capability of driving the output y(t) to some desired 
constant value y d on command. Generally, this can only be accomplished by driving 
the control vector u(t) to some constant value Up which is required to satisfy the 
steady-state solution of (1.1-2), i.e., 

FXy + Gu r = 0 (l.3) 

y d = HXy (l.H) 

where the equilibrium point represented by the pair [x r ,u r ] is unknown. A control law 
that will simultaneously drive -the output y(t) toward. y d and the control u(t) toward 
u,- for an arbitrary choice of both y d and u r may be obtained by minimizing a scalar 
cost functional of the form 

V[x(t 0 ) ,t 0 ,T;u( • ) ] = [ II y ( t ) - y d li* + Mt) - V^jdt (1.5) 

t o 

where the matrices Q and R are symmetric, positive definite, constant, and of appro- 
priate dimension. The symbol llyll^ represents y’Qy where the prime indicates 
transpose. 

The cost functional (1.5) contains an unknown constant vector u r . In addition to the 
control law that minimizes (1.5), the vector Uj. must be determined so that the opti- 
mal scalar cost is finite as the final time T approaches infinity. This will ensure 
that the resulting optimal control law will drive the output y(t) to y d . 


This tracking problem can be put in the form of an output regulator problem by defining 
new variables: 


and new matrices : 


Then for a plant represented by 

x(t) = Fx(t ) + Gu(t ) ; x'(t 0 ) = [x Q » Uj. y d ] (l.9) 

y(t) = Hx(t ) (1.10) 



u(t) 

= u(t) 

- % 

9 

y(t) = y(t) - y d 

(1.6) 



x'(t) 

= [x(t) 

^ y d ] 

(1.7) 

[F 

G 

°1 


"g‘ 



0 

0 

0 ; 

G = 

0 

; H = [H 0 -I] . 

(1.8) 

Lo 

0 

oj 


. 0 . 
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( 1 . 11 ) 


the control law that minimizes a quadratic cost functional 

V(x(t Q ) ,t 0 ,T;u( • ) ) = [ly(t)!^+ l|u(t)fl^]dt 

is given by ° 


u#(t) = -R _1 G'P(t,T)x(t) 


where P(t,T) is the solution of the matrix Riccati equation 

p(t) = -P(t )F - F'P(t) + P(t)GR" 1 G'P(t) - H'QH ; P(T) = 0 . 


( 1 . 12 ) 

(1.13) 


The optimal value of the scalar cost functional (l.ll) using the optimal control law 
(1.12) is given by 

v*(x(t 0 ),t 0 ,T) = | ilx(t 0 )!!p (to}T) . (1M 

For all T, the solution P(t,T) exists, is unique, is nonnegative definite symmetric, 
and is defined for all t < T.( 3 ) 


If P(t,T) is partitioned as 

■p(t,T) Pi 2 (t,T) Pi 3 (t,T) 

P(t,T) = P 12 (t ,T) P 22 (t,T) P 23 (t,T) 

-Pi 3 (t ,T) P 23 (t ,T) P 33 (t,T)J 

the optimal control law (1.12) can be written as 

u*(t) = -R -1 G'P(t ,T)x(t ) - R- 1 G*[P 12 (t,T)u r + Pi 3 (t,T)y d ] . 
The matrix Riccati equation (1.13) becomes 


(1.15) 


(1.16) 


P(t) = -P(t)F - F'P(t) + P(t)GR _1 G'P(t) - H'QH ; P(T) = 0 (l.lT) 

P i2 (t) = -(F - GR -1 G'P(t ) ) 'Pi 2 (t ) - P(t)G ; Pi 2 (T) = 0 (l.l8) 

P i3 (t) = -(F - GR“ 1 G'P(t)) , P 13 (t) + H’Q ; Pi 3 (T) = 0 (1.19) 

P 22 (t) = P 1 2 ( t ) GR~ 1 G'P 12 (t) - Pi 2 (t)G - G'Pi 2 (t) ; P 22 (T) = 0 (l.20) 

P 23 (t) = -G'P 13 (t) + Pi 2 (t)GR- 1 G , Pi 3 (t) ; P 23 (T) = 0 (l.2l) 

P 33 (t) = Pi 3 (t)GR- 1 G’P 13 (t) - Q ; P 33 (T) = 0 (1.22) 

and the optimal scalar - cost (l.lU) becomes 

V*(x(t 0 ),t 0 ,T) = |j + 2x^[P 12 (t 0 ,T)u r + Pi 3 (t 0 ,T)y d ] + u;P 22 (t D ,T)u. 

+u'P 23 (t 0 ,T)y d + y^P 23 (t Q ,T)u r + y«P 33 (t 0 ,T)y d | . (l.23) 

If we define 

b(t) = -P 12 (t)u r - Pi 3 (t)y d (1.2U) 

c(t) = u^P 22 (t)u r + u^P 23 (t)y d + y^P 23 (t)u r + y d P 33 (t)y d (l.25) 

equations (1.18-19) can be combined to form 

B(t) = -(F - GR _1 G'P(t) ) ’b(t) - (H’Qy d - P(t)Gu r ) ; b(T) = 0 (1.26) 

and (1.20-22) can be combined to form 

c(t ) = b'(t)GR” 1 G , b(t) - y *Qy + b'(t)Gu + u'G'b(t) ; c(T) = 0 . (1.27) 

d d r r 

Thus, solution of the matrix Riccati equation (1.13) is equivalent to simultaneous solu- 
tion of (1.17), (1.26), and (1.27). Since u#(t) = u*(t) - Uj., the optimal control law 
for the plant (1.1-2) is given by 

u *(t) = -R _1 G'P(t,T)x(t) + R" 1 G'b(t,T) + u r (1.28) 
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where b(t,T) is the solution of ( 1 . 26 ). The optimal scalar cost using (1.28) is 
given by ~ 

V*(x(t 0 ),t 0 ,T) = | [llx 0 llp (to>T) - 2xJb(t Q ,T) + c(t 0 ,T)] (1.29) 

where c(t,T) is the solution of (1.27). Equation (1.17) is the matrix Riccati equa- 
tion for the output regulator problem for the plant (1.1-2). 

The opti ma l control law (1.28) and the resulting optimal scalar cost (1.29) are func- 
- tions of the final time T and the arbitrary vector- Uj.. . The limiting solution as 
T -*• * can be shown to have the following properties: 

P = lim P(t,T) = lim P(t,T) (1.30) 

T-x® t - * - — 00 

exists, is constant, and can be found by obtaining the steady-state solution of (1.17). 

b = lim b(t,T) = lim b(t,T) (l.3l) 

T-Xxj t->-— 00 

exists, is constant, and can be found by obtaining the steady-state solution of the 
linear, first-order differential equation 

b(t) = (F - GR" 1 G»P)'b(t) + (H’Qy d - PGu^) ; b(t Q ) = 0 . (1.32) 

For arbitrary ul, , solution of (l.32) yields 

b = Px (1.33) 

r 

such that satisfies 

r 

P(Fx + Gu ) + H*Q(Hx - y ) = 0 . (1.3U) 

r r r a 

A bounded optimal scalar cost (1.29) as T -*■ <=° requires that 6(t) = 0 or c(t) = c a 
constant, yielding the optimality condition 

b , GR“ 1 G'b - y'Qy. + b'Gu + u'G'b = 0 (1.35) 

d d r r 

i.e., if u r satisfies (l.35) where b is obtained from (l.32), then the control law 

u*(t) = -R -1 G'Px(t) + R -1 G’b + Up (1.36) 

is optimal, the controlled plant 

x(t) = (F - GR _1 G'P)x(t) + Gu* (1.37) 

u* = R _1 G'b + Uy (1.38) 

is asymptotically stable for all t, and the optimal feedforward control law u* will 
drive the output y(t) to y d . The vectors x r , Uj., and y d satisfy (l.3-^) where 
x r = P -1 b, the constant c is given by 

c = x’Px (l.39) 

r r 

and the optimal scalar cost as T -* 00 is given by 

V*(x(t Q ),t 0 ,») = | [||x 0 ijp - 2x^b + c] = \ ||x 0 - xjjp . (l.J+O) 

The problem now is: given a desired constant output vector y d , find the corresponding 

constant required control vector Uy that satisfies (l.35) so that the optimal feed- 
forward control law (1.38) can be computed. The problem has been solved by developing 
the iterative algorithm given in the following steps. 

A. Solve (1.17) for the steady-state value P. 

B. Set ul. = 0. 

C. Solve (1.32) for the steady-state value b. 

D. Substitute Uy and b into (1.37-38) and solve for the steady-state value x. 

E. Substitute x into (1.2) and solve for y. 
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F. Compare this value of — y with y^. If a significant difference exists, compute u* 
hy substituting x, b, and u r into (1.36), set Up = u*, and return to step £7. 

The iteration is terminated when step F shows y to be sufficiently close to y^. The 
feedforward control law (1.38) used to achieve this result is considered optimal. 

A sufficient condition for convergence of the algorithm is that the magnitude of the 
eigenvalues of the matrix I - HrH'Q are less than unity where 

T = E’PGR^G’PE (l.4l) 

and 

E = -[P(F - GR“ X G'P) ] _1 . (1.42) 

2. SOLUTION OF THE OBSERVER EQUATION: TF - AT = B 

Let r components of the n-dimensional state vector x(t) of the plant represented 
by (l.l) be considered unmeasurable, where 

1 S r js n-1 . (2.1) 

Let z(t) be an r-dimensional column vector defined as the output of a filter of the 
form 

z(t) = Az(t) + Bx(t) + Cu(t) ; z(t Q ) = z 0 • (2.2) 

The matrices A, B, and C are constant, rectangular, and of appropriate dimension. 

It has been shown( 4 ) that if a constant, rectangular transformation matrix T can be 
found satisfying the matrix equation 

TF - AT = B (2.3) 

where A and F have no common eigenvalues, then, after a time associated with the 
eigenvalues of the matrix A, the state x(t) of the plant (l.l) is related to the 
output z(t) of the filter (2.2) by the transformation 

z(t) = Tx(t) (2.4) 

and, in addition, the matrix C in (2.2) is given by 

C = TG . (2.5) 

Therefore, if we could determine the matrices T, A, and B that satisfy (2.3), we 
could solve (2.4) for an estimate of those states which sire considered unmeasurable in 
terms of the measurable states and the filter output. We could then construct a state 
vector from the measurable states and the estimates of the unmeasurable states which 
could be used in the feedback control law for the output regulator. 

The observer filter (2.2) cannot use any of the unmeasurable states as input. There- 
fore, the columns of the matrix B associated with the unmeasurable states must have 
all zero elements. It would also be desirable if the solutions of (2.3) were such that 
the eigenvalues of A were large and negative to ensure quick response, and that the 
absolute value of the elements of T, B, and C were small enough to be practical. 


Since we know nothing about the elements of T, A, and B, except that certain columns 
of B must be zero, there are many more unknowns than equations. Certain constraints 
are therefore imposed to facilitate the computational task. 


Restrict the matrix A to be diagonal. This restriction greatly simplifies the filter 
equations for multidimensional filters by eliminating cross-coupling between the ele- 
ments of the filter vector. Temporarily assume that all eigenvalues of the matrix A 
are equal and denoted by a. Then 

A = al (2.6) 


and (2.3) becomes 


T[al - F] + B = 0 . 


(2.7) 
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The equation for the 


where b, and c. are 


i 


i 


jth_ element of the filter becomes 
z^(t) = aZj(t) + bjx(t) + c^u(t) 
the jth rows of B and C, respectively. 


( 2 . 8 ) 


The problem therefore has been reduced to that of solving a matrix equation of the 
form (2.7) to obtain r simple, uncoupled, first-order filters of the form (2.8). 


Assume that the kth element of the state vector x(t) is to be estimated primarily 
by z j(t). We would thus like t^ (the kth component of tj, which is in turn the 
jth row of T) to take on the value unity, and all other components of tj and all 
components of bj and Cj to be small in magnitude. 


To achieve a -unique solution of (2.7), r cost functionals of the following form will 
be minimized subject to the constraint of (2.7): 


WV = 2 


[(t j - V (t j - S’ 


( V E ° (b j )] 


(2.9) 


where R 0 is a constant diagonal matrix and 
nent which is unity. 



is zero except for the kth compo- 


A standard LaGrange multiplier minimization yields a solution tj and b j . To assure 
the desired zero columns in B, it will be assumed that the inverse R^ 1 has zeros at 
the appropriate diagonal components; i.e., R 0 itself assigns infinite weight to the 
columns that must be zero. After the minimization procedure, the rows t j , b j , and Cj 
are normalized to assure that t^v is unity. The normalized matrices are denoted by 
(~). J 


The preceding is implemented with the following algorithm. The algorithm incorporates 
the proper choice of the eigenvalues of the matrix A and the construction of all r 
filters simultaneously. 


A. Given a plant represented by (l.l), select r elements of the state vector which 
are to be assumed unmeasurable, where 

1 < r < n-1 . 

B. Set the corresponding elements of the diagonal matrix Rq 1 to zero and the remain- 
ing elements to unity. 

C. Construct a projection matrix T 0 such that if z(t) were defined as 

z ( t ) = T 0 x(t) (2.10) 

then the elements of the vector z(t) would be the unmeasurable states. (T Q has one 
unity element in each row and zeros otherwise.) 

D. Select a positive gain limit such that if the absolute value of any element of the 
solution is greater than this gain limit the solution is to be considered impractical. 

E. Choose the parameter a to be some small negative value whose magnitude is approxi- 
mately equal to the magnitude of the smallest eigenvalue of the matrix F. 

F. Construct the matrix 

M = [al - F] (2.11) 


G. 


Carry out the minimization. 


Construct and solve the following matrix equations: 



( 2 . 12 ) 


B' = 

C = TG 
A = la 
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(2.13) 

(2.14) 

(2.15) 



which yields the following observer 

z(t) = Tx(t) (2.l6) 

i(t ) = Az(t) + Bx(t ) + Cu(t) (2.17) 

H. Normalize the observer. For each row of T os find the unity element and then select 
the corresponding element of the same row of T (i.e., tj^K Multiply this same row of 
T, B, and C by the inverse of this element (i.e., 1/t , yielding the normalized 

(2.18) 

(2.19) 


observer 

z n (t) = Tx(t ) 

z n (t) = Az n (t) + 5x(t) + Cu(t) 

where z n (t) is used for the normalized variable. The matrices T, B, and C 

TF - AT = B 
C = TG 


and 


satisfy 

( 2 . 20 ) 

( 2 . 21 ) 


X. Consider the rows of the matrices A, T, B, and C to be solutions. Compare the 
absolute value of the elements of each solution with the gain limit. If the absolute 
value of each element of a solution is less than or equal to the gain limit, store 
that solution. If the absolute value of any element of a solution is greater than the 
gain limit, do not store that solution. If the absolute value of at least one element 
in all solutions is greater than the gain limit, print the stored solutions and stop. 

If not, increment the parameter a to some larger negative value and return to step F. 

If the matrix inverse indicated in (2.12) exists, each iteration will produce r 
independent solutions as defined in step J of the algorithm. ’ The iterative procedure 
will cause the r solutions with the quickest response time, whose elements do not 
exceed the gain limit in magnitude, to be stored. All other solutions are discarded. 
The diagonal elements of the resulting stored matrix A are not generally equal. 

An illustrative example of the use of the algorithms will now be presented. 

3. ILLUSTRATIVE EXAMPLE 

The aircraft chosen to illustrate the procedure presented here is the Augmentor 
Wing Jet STOL Research Aircraft. For longitudinal flight, in addition to the 
usual elevator and throttle, this aircraft has vectored thrust provided by four Pegasus 
nozzles that direct the hot gas exhaust from the engines. It also has augmentor blow- 
ing ejector flaps through which-the cold airflow from the engines is ducted. 

For a 60-knot approach on a -7.5° flight-path angle, linearized longitudinal perturba- 
tion equations were obtained of the form (1.1-2) where 

x'(t) = [v(t) a(t ) 6(t) q(t ) ] 


u’(t) = [S E (t) 
y * (t ) = [v R (t) 


6 N (t) 


« T (t) 


yt» 


a(t) e(t) y(t)] 


(3.1) 

(3.2) 

(3.3) 


and 

v(t) = variation in the forward component of the airspeed, ft/sec 
a(t) = variation in angle of attack, deg 
6(t) = variation in pitch attitude, deg 
q(t) = variation in pitch rate, deg/sec 
5g(t) = variation in elevator angle, deg 
6jj(t) = variation in nozzle angle, deg 
<S,p(t ) = variation in throttle angle, deg 
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6p(t) = variation in flap angle, deg 

v K (t) = variation in the forward component of the airspeed, knots 
y(t) = variation in the flight-path angle, deg 


The matrices 


F, G, and H are 


[F] 


-0.05479 

-0.1453 

0 

_ 0.08273 


0.2606 

-0.5078 

0 

-0.07091 


-0.5598 

0.02461 

0 

-0.00965 


-0.0880" 

0.9521 

1.0 

-1.352 . 


[G] 


" 0. 004216 
-0.04573 
0 

-1.302 


-0.08378 

0 

0 

- 0.08190 


-0.01867 

- 0.5028 

0 

0.2676 


-0.1031" 

-0.09994 

0 

0.2167. 


[H] 


0.5925 

0 

0 

.0 


0 0 0 " 

1.0 0 0 

0 1.0 0 * 

- 1.0 1.0 0 . 


(3.4) 


(3.5) 


(3.6) 


The matrices Q and R" 1 were chosen to be 



‘0.4 

0 

0 

0 " 


“0.625 

0 

0 

o' 

[Ql - 

0 

0 

1.6 

0 

0 

1.6 

0 

0 

; [jrM - 

0 

0 

10 

0 

0 

2.5 

0 

0 


_0 

0 

0 

10. 


.0 

0 

0 

0. 


(3.7) 


For computational convenience R -1 was taken as a 4 x 4 matrix with a zero row and a 
zero column since the flap was not used for feedback control. The matrix Riccati equa- 
tion was solved. The control law u(t) » -R” 1 G'Px(t) where 


[r” 1 g , p] 


0.1009 

-0.6807 

0.07422 
. 0 


-0.2650 

-0.03284 

-4.392 

0 


-0.9779 

- 1.049 

4.953 

0 


- 0.6761 

-0.5776 

0.3067 

0 


(3.8) 


gave reasonably good response to initial condition errors and so was taken as an accep- 
table basis for the design of optimal feedforward control laws. 


A number of optimal feedforward- control laws were computed for various desired steady- 
state outputs y d . The desired output criterion for each case is shown in table 1. 


Table 1.- Steady-State Output Design Criteria 


s sQase 

y<i\ 

■ 

2 

3 

■ 

5 

6 

■ 

8 

9 

10 

11 

V v 


0 


0 

Free 

Free 

Free 

Free 

10.0 


10.0 


0 

-1.0 


Free 

0 

-1.0 

Free 

Free 

0 

9 

Free 

a 

0* 

1.0 

0 

R989 

1.0 

1.0 

0 

Free 

1.0 

0 


Free 

a 

Yd 

1.0 

Ba 


Free 

1.0 

1.0 

1.0 

Free 

0 

0 

Free 


Note that, in some cases, some of the elements of the output vector were allowed to 
remain free while others were to be driven to a specific value. This was achieved by 
starting with the free elements of y d equal to zero in step B of the algorithm and 
by replacing the free elements of y d with the corresponding elements of y in step 

F. 
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The desired output vector y^ and the required control vector Uj. were computed for 
each case. The resulting- steady-state output and control vector pairs for each case 
are shown in tables 2 and 3, respectively. The corresponding optimal feedforward con- 
trol law u* for each case is shown in table 4. 
c 

Table 2.- Desired Steady-State Output Vectors 


*N^ase 

1 

2 

3 

b 

5 

6 

7 

8 

9 

10 

11 

v k 

0 

0 

0 

0 

enroii 

-.3120 

-.5275 

-.8907 

10.0 

10.0 

10.0 

K d 

0 

-1.0 

-.8095 

.9332 

■SIS 

-1.0 

-.7513 

.9210 

0 

-1.103 

-.9281* 

6 d 

1.0 

0 

.1905 

1.0 

• 

0 

.2487 

1.0 

0 

-1.103 

-1.144 

^d 

1.0 

1.0 

1.0 

.06676 

ESI 

1.0 

1.0 

.07916 

0 

0 

-.2157 


Table 3.- Required Steady-State Control Vectors 


XCase 

V\ 

1 

2 

3 

m 

5 

6 

■ 

8 

9 

10 

11 

6 E r 

6 N r 

6 T r 

6 F r 

o.4i4 

-6.663 

.011 

0 

0.U61 

-3.303 

.968 

0 

0.1*52 

-3.943 

.786 

0 

-.016 

-3.581 

-.892 

0 

.334 

-5.497 

.593 

0 

.44o 

-2.994 

1.122 

0 

.413 

-3.6l6 

.991 

0 

-.071 

-2.740 

-.41*0 

0 

.680 

-9.907 

-4.941 

0 

.732 

-6.201 

-3.886 

0 

.635 

-5.351 

-4.056 

0 


Table 4.- Optimal Feedforward Control Laws 


ng 

1 

2 

3 

n 

5 

6 

7 

8 

9 

10 

11 

a* 

°E C 

-.564 

.726 

.480 

-1.242 

-.845 

.652 

I 

-1.445 

2.38 

3.80 

3.70 

«I C 

-7.712 

-3.270 

-4.117 

-4.660 

-5.191 

-2.603 


-2.796 

-21.40 

-16.50 

-15.61 

r* 

6 t c 

4.964 

5.360 

5.285 

-.038 

5.399 

5.475 

5.456 

.357 

-3.68 

-3.25 

-4.39 

t* 

°F C 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


For each case, the optimal feedforward commands (table 4) will drive the output y(t) 
qP the controlled plant (l.3T) to the desired steady— state output (table 2) with the 
resulting steady-state control deflections (table 3). Transient responses for case 9, 
case 2, and case 1 are shown in figures 1, 2, and 3, respectively, where the elevator, 
nozzle, and throttle were assumed to have first-order actuator dynamics with time con- 
stants of 0.1, 0.0667, and 2.5 seconds, respectively. These figures show that the 
steady— state responses are exactly as desired, with quite reasonable dynamics as well. 


Let us assume that v(t) and a(t) are to be considered unmeasurable. Then if we 
construct 



0 0 

0 0 

1.0 0 

0 1 . 0 . 



0 

1.0 


0 0*1 

0 oj 


(3.9) 


and use a gain limit of 10, an 


observer of the form (2.18-19) 


is obtained where 
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U:°o44 i:® 53 

-0.7202 

3 . l4ll 

(3.10) 

0.1679 

6.865J 

w = [1' 5 h] 



(3.11) 

r ~! [o o -0.8308 

LBJ " Lo 0 0.3695 

1.136 I 

-0.5744J 


(3.12) 

r ~, f-4.308 -0.3411 

LLJ “ L-8.986 -0.5284 

-I.618 

0.092641 

(3.13) 

1.342 

1.429 J 


Figure 4 shows the transient response for case 1 using this observer as dynamic compen- 
sation, and using only pitch attitude and pitch rate for feedback. The same actuator 
dynamics were assumed. Comparison of figure 4 with figure 3 shows an almost identical 
response. 

4. CONCLUSIONS 

Simple algorithms for the construction of optimal feedforward control laws and optimal 
observers have been presented. These algorithms can be used to obtain control laws 
that behave as stable regulators for the. plant in the absence of command inputs, that 
will drive the output of the plant to a desired constant value in response to a step 
command input, and that require only selected elements of the state vector for feedback. 
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Figure 1.- Airspeed command input with entire state measurable. 
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Figure 2.- Direct lift command input with entire state measurable. 
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Figure 3.- Flight-path command input with entire state measurable. 
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Figure 4.- Flight-path command input with second-order observer. 
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